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Abstract: We present the leading experimental constraints on snpersymmetric models 
with R-parity violation (RPV) and a long-lived lightest snperpartner (LSP). We consider 
both the well-motivated dynamical RPV scenario as well as the conventional holomorphic 
RPV operators. Gnided by natnralness, we stndy the cases of stop, glnino, and higgsino 
LSPs with several possible leading decay channels in each case. The CMS displaced dijet 
and the ATLAS mnltitrack displaced vertex searches have been fnlly recast, with all cnts 
and vertex reconstrnction algorithms applied. Heavy charged stable particle searches by 
CMS are also applied. In addition, we consider representative bonnds for prompt LSP 
decays that are directly applicable. Onr main resnlts are exclnsion plots in the mnsp — tlsp 
plane for the varions scenarios. We find that the natnral parameter space {mi < 800 GeV, 
mg < 1500 GeV, mg < 800 GeV) is exclnded for a long-lived LSP (tlsp ^ 1 mm). 
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1 Introduction 

For the past three decades, supersymmetry has been widely considered the top candidate 
for completing the Standard Model (SM) of particle physics. The expectation was that if 
supersymmetry provides the solution to the naturalness problem, the superpartners should 
lie somewhat below the TeV-scale, the energy regime probed by the Large Hadron Collider 
(LHC). However, Run-I of the LHC did not yield any evidence for supersymmetry [1-5], 
challenging our preconception of naturalness. 

LHC searches for supersymmetry typically imply the superpartner masses in the mini¬ 
mal supersymmetric extension of the Standard Model (MSSM) to be above 1 TeV, rendering 
most of the parameter space unnatural. Thus if supersymmetry is to remain natural, it 
must manifest itself non-minimally. This has led to many new ideas of non-minimal models 
and renewed interest in older ideas, such as R-parity violation (RPV) [6-12]. Many of the 
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constraints from standard superpartner searches can be evaded if the assumption of exact 
R-parity conservation is relaxed (see e.g. [12-25] and references therein). In particular, 
if the LSP decays inside the detector, then supersymmetry searches which require large 
missing transverse energy (MET) will not be applicable to these models. 

At first sight the inclusion of RPV operators appears to lead to a new set of problems: 
each of the plethora of new couplings introduced has to remain small in order to suppress 
baryon- and lepton-violating processes. This suggests that there should be a systematic 
suppression of all RPV operators. A simple dynamical reason behind such a suppression 
would be if the visible sector (comprised of the SM fields and their superpartners) is R- 
parity conserving, while R-parity is violated in a hidden sector. The mediation of R-parity 
violation via heavy vector-like messengers will give rise to effective RPV operators in the 
visible sector, suppressed by the mass scale of the heavy messengers and possibly also by 
light fermion masses or the supersymmetry breaking scale. Within this general framework, 
called dynamical R-parity violation (dRPV) [26, 27], all R-parity violation in the visible 
sector originates from effective higher dimensional operators. 

Since all RPV operators are dynamically suppressed, one needs to identify the leading 
RPV operators in the visible sector. As shown in [26, 27], this identification depends on the 
details of the mediating dynamics and, in particular, on the U{\)b-l and U{1)r charges 
of the fields responsible for the breaking of RPV. One finds that while in some cases the 
standard holomorphic superpotential RPV operators, 

IA^RPV — T T T ^ijk^i^j^ki (^'l) 

may dominate, typically the leading R-parity-violating operators will originate from the 
non-holomorphic, non-renormalizable Kahler terms, 

^nhRPV = ^ ■ (1-^) 

The interactions from (1.2) are suppressed either by the SUSY-breaking scale or by the 
fermion masses. The structure of the novel R-parity-violating operators (1.2) and other 
non-standard RPV operators lead to new and distinct collider signatures. 

Both the standard RPV couplings (1.1) and the non-holomorphic couplings (1.2) must 
be small. This often leads to LSP lifetimes of ctlsp ^ mm, corresponding to particles 
which are long-lived on collider scales. Thus the LSP can appear in the detector as an 
exotic non-prompt decay without necessarily producing missing energy. Both ATLAS and 
CMS have performed searches for such non-prompt decays, including searches for displaced 
vertices, decaying R-hadrons, and heavy stable charged particles (HSCP). The aim of this 
paper is to explore the final states allowed by the operators of (1.1) and (1.2) and establish 
the resulting bounds in the LSP mass vs. lifetime parameter space. Our focus is on the 
cases motivated by naturalness, where the LSP is a stop squark, a gluino or a Higgsino. We 
recast the main ATLAS [28] and CMS [29] searches for displaced vertices for the various 
R-parity violating scenarios. Some of the selected final state signatures are unique to the 
non-holomorphic operators (1.2) {e.g. t —?■ 66), while others have bounds which apply 
to both holomorphic and non-holomorphic RPV operators. HSCP searches by CMS are 
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reinterpreted for the scenario where the LSP is unstable. We also present some results for 
prompt searches by ATLAS and CMS to determine whether there exists unexplored gaps 
in the parameter space for short-lifetime LSPs. 

While this work was nearing completion, a closely related analysis [30] by Zhen Liu and 
Brock Tweedie appeared which has a significant overlap with our paper. While our work 
studies the various RPV topologies, [30] focuses on various supersymmetric scenarios. Our 
work also incorporates the recasting of the most recent ATLAS displaced vertex search [28] 
and uses slightly different tracking efficiencies and vertex reconstruction procedures. 

The paper is structured as follows. In Sec. 2 we present the models and parameter 
space which we consider for our analysis. In Sec. 3 we briefly describe the experimental 
searches considered in our analysis. In Sec. 4 we present our results in the form of exclusion 
plots in the LSP mass vs. lifetime parameter space. We conclude in Sec. 5. In Apps. A and 
B we detail the analysis techniques implemented in the various experimental searches we 
have recast and present our procedure used for each channel. 

2 Theory and Models 

The aim of this project is to explore the RPV scenarios that are most relevant for maintain¬ 
ing a natural Higgs sector within the supersymmetric theory. This requires the stop be well 
below the TeV scale, the gluino lie not too far above the stop, and the Higgsinos to be light 
as well [31]. Thus we study models with light stops, gluinos, or Higgsinos below the TeV 
scale. The LHC phenomenology will strongly depend on the identity of the LSP and the 
nature of the RPV operator responsible for its decay. In accordance with the naturalness 
expectations, we consider three types of LSP particles: stops (t), gluinos {g), and Higginos 
(H). 

Since a main focus of this paper is the study of the non-holomorphic RPV operators 
of the dRPV scenario in (1.2), we concentrate on decays that predict final outgoing jets (in 
addition to other possible final-state particles). The reason for this is that the SM gauge 
structure of the non-holomorphic operators (1.2) differs from the holomorphic operators 
(1.1). In particular there is no non-holomorphic cubic term that contains more than one 
lepton (there is no analog of the lie coupling), thus a jet is expected in every final state 
involving these operators. An exception which we do not consider could be the case for a 
Higgsino LSP if the breaking of R-parity occurs via bi-linear terms or via the holomorphic 
lie operator. 

For the various LSPs, the decay channels considered in this paper are summarized 
in Table 1. We only consider direct production of LSPs. We study three possible decay 
channels for the stop: 

• t ^ dd! where d, d! are down-type SM quarks. Such decays arise either from the 
\" holomorphic RPV operator in Eq. (1.1) or from the rj” operator of the non- 
holomorphic dRPV operators in Eq. (1.2). While decays to only third generation 
particles t —)• 66 are not allowed with the holomorphic X” operator due to the an¬ 
tisymmetry of the color contraction, such decays are allowed and are expected to 
dominate in the case of the rj" operator in Eq. (1.2) . 


- 3 - 


Topologies 

LSP 

Decay 

Operator 

Results 


dd' 

A", g" 

Fig. 1 

i 

U V 

g' 

Fig. 2 


dl+ 

A', g 

Fig. 3 


tdd' + c.c 

A", g" 

Fig. 4 

9 

tuv + c.c 

g' 

Fig. 4 


td£~ + c.c 

A', g 

Fig. 5 


(t/b) dd' + c.c 

A", g" 

Fig. 6 


(t/b) uv + c.c 

g' 

Fig. 6 


{t/b) di~ + c.c 

A', g 

Fig. 6 


Table 1. Summary of the various LSPs and their decay channels considered in this paper. The 
third column denotes the RPV operators from (1.1) or (1.2) that can give rise to the decays and 
the final column points to the relevant exclnsion plot in Section 4. 

• t —)• nP where u is an up-type quark and P is an anti-neutrino. This channel occurs 
for the r]' operator of Eq. (1.2) but is suppressed by the neutrino mass. However, this 
decay can be induced by effective Kahler potential terms of the form JC D (V^qu^ 
which were shown in [27] to be generated in a simple model of dRPV. The decay 
can also be induced via superpotential terms of the form W D ihuhuQU and via 
neutrino-higgsino mixing. 

• t —7- dl^ where is an anti-lepton. The responsible operator may be of Eq. (1.1) 
or r] of Eq. (1.2). 

The chiral suppression of the decays from the non-holomorphic operators and the typ¬ 
ical flavor structure of RPV operators strongly suggests that stop decays to heavy quarks 
dominates over decays to light quarks. Thus we expect that third generation hnal states 
will dominate both for the case of stop ESP and for the other possibilities discussed below. 
Since the stop is assumed to be the lightest squark, the gluino decays via an off-shell stop 
and the hnal states will correspond to the above listed hnal states for stop decay with an 
additional top quark, g —>■ tdd' + c.c., g —)> tuv + c.c. and g —>■ tdl~ + c.c., with third gen¬ 
eration hnal states again expected to dominate. Similarly, a neutral Higgsino LSP decays 
via an off-shell stop and its topologies are identical to those of the gluino LSP decays. A 
charged Higgsino will instead have a bottom quark replacing the top quark in the respective 
hnal states. 

3 Experimental Searches 

In this section, we give a brief overview of the searches we have used for our study. A more 
detailed discussion of the individual searches, as well as our methods of applying them to 
RPV models, can be found in Appendices A and B. 
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The LSP may decay promptly at the primary vertex (PV) or may traverse a finite 
distance before decaying, depending on the size of the corresponding RPV coupling. While 
we mainly focus on displaced LSP decays due to RPV interactions, we consider the limits 
obtained from prompt and HSCP studies as well (corresponding to very short or long LSP 
lifetimes respectively). 

3.1 Searches for Displaced Vertices 

LSPs which decay in the inner detector can be reconstructed as displaced vertices by both 
the ATLAS and CMS trackers. The SM background for these searches is negligible. Long 
lifetimes, for which only a few of the produced LSPs decay within the tracker, can be ex¬ 
cluded, ranging from 1 to lO”^ mm proper lifetimes. The searches presented in this paper 
look for displaced vertices associated with one of several signatures used for triggering, 
including jets, leptons, and MET. We recast displaced searches by both ATLAS [28] and 
CMS [29], which together place strong constraints on all of the final states considered in this 
paper. The full detail of the simulation and procedure, including cuts, tracking efficiencies, 
and vertex reconstruction, are given in Appendices A and B. 

ATLAS DV+/i/e/jets/MET. 

The ATLAS search for long lived particles at \/8 TeV [28] considers four multitrack signa¬ 
tures with a luminosity of 20.3 fb~^ in which a displaced vertex (DV) with at least 5 tracks 
is accompanied by one of the following: (1) a high px muon (2) a high pT electron (3) jets 
or (4) missing transverse energy (MET). These searches are background free. The muon 
search has been improved from the previous study at 8 TeV [32] by increasing the detector 
volume used for reconstructing vertices. 

CMS Displaced Dijet. 

The CMS collaboration has performed a search for dijets originating from displaced vertices 
at ^/s = 8 TeV with 18.6 fb“^ of data [29]. The analysis requires two hard jets consistent 
with a dijet produced at a secondary displaced vertex. The constraints from this search are 
not limited to models with dijets. Isolated leptons are treated as jets, and three-jet final 
states have significant efficiency to be reconstructed as dijets, giving sensitivity to a variety 
of models. 

3.2 Searches for Stable Particles 

Long-lived colored particles, such as squarks or gluinos, form composite objects with SM 
particles, called R-hadrons. For long enough lifetimes, the LSP may hadronize and traverse 
the whole detector before decaying. Searches for slowly moving charged particles constrain 
the production rate of a stable LSP. These searches also place strong limits on unstable 
LSPs with lifetimes long enough for a fraction of the particles to traverse the detector before 
decaying. 
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CMS Heavy Stable Charged Particle Search. 

CMS has performed a search for heavy stable charged particles (HSCP) at ^/s = 8 TeV with 
an integrated luminosity of 18.8 fb“^ [33]. HSCPs have velocities much less than the speed 
of light leading to misidentihcation in the particle flow reconstruction algorithm used by 
CMS which assumes SM particles with r; ~ c. Additionally, tracks associated to the HSCP 
typically have large energy loss {dE/dx) due to interactions with the tracker. Thus, HSCPs 
are identihed by a longer time-of-flight to the muon system or by tracks with anomalous 
energy loss. i?-hadrons can be either charged or neutral, but interactions with the detector 
material may "charge-strip" the charged i?-hadrons, converting them into neutral hadrons 
by the time they reach the outer muon system. To account for this possibility, CMS 
considers two models of the R-hadrons’ strong interactions with the material in the detector. 
The hrst, referred to as the cloud-model, assumes that the HSCP is surrounded by a cloud 
of light SM colored particles. The second is a charge-suppressed model, in which charged 
i?-hadrons quickly become neutral by interacting with the detector material and are often 
neutral before reaching the muon system. 

3.3 Prompt Searches 

These searches are relevant for LSPs with small proper lifetimes, taken conservatively to be 
below 1 mm in this study, for which the decay products of the LSP point back to the primary 
production vertex. We have not recast prompt searches. Instead, we have used existing 
bounds which are directly applicable to the particular scenarios under consideration. The 
searches used include the paired dijet resonance search at CMS [34] for the fully hadronic 
prompt stop decays, leptoquark searches [35-43] from the Tevatron and the LHC for t —>■ dl~^ 
decays, R-parity conserving supersymmetric searches [1-5, 44] for t —>■ tN, c —>■ cN, g —)> tiN 
(for vanishing neutralino mass) for decays involving neutrinos, and the R-parity violating 
ATLAS gluino search [45]. 

4 Results 

In this section we present our results for the bounds on the LSP masses and lifetimes in 
the various scenarios considered and summarised in Table 1. Our results follow from a 
detailed MC study which includes all channels and the application of the cuts and the 
reconstruction procedures, designed to correspond to the various experimental displaced 
vertex searches. We postpone the details and discussion of these procedures to Appendix A. 
The HSCP search was recast at the parton level in order to hnd the efficiency for an unstable 
LSP to traverse the detector before decaying. Supersymmetric production cross sections 
were obtained from the LHC SUSY cross section working group [46, 47] when applicable. 
Higgsino direct production was calculated at tree-level in Madgraph, and an NLO k-factor of 
1.6 was applied for all masses [48, 49]. Exclusion bounds from prompt searches are directly 
taken from the experimental results, and were not recast, implying that only a small fraction 
of the searches relevant to prompt decays have been included. To obtain further bounds in 
the prompt case requires recasting those searches and is beyond the scope of this paper. 
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Figure 1. 95% CL exclusion curves for fully hadronic stop decays to two down-type quarks. In 
the upper left plot we show i —>■ dd, in the upper right plot i —)■ db, and in the bottom plot 
i ^ bb is presented. While all displaced searches of ATLAS [28] and CMS [29] have been recast, 
we only display the curves yielding the strongest constraints. Stable searches for HSCPs [33] were 
recast for an unstable R-hadron. The prompt searches here correspond to the CMS dijet resonance 
search [34] . 


4.1 Stop LSP 

For the case of a stop LSP, we have considered the three final state topologies, as given in 
Table 1 . Only direct production of the stop has been simulated and the gluino was assumed 
to be decoupled for the determination of production cross section. Despite the change in 
kinematics, we expect the gluino production to mostly affect the overall stop production 
rate, and therefore the limits presented here can be used to estimate the limits for the 
non-decoupled gluino case. 

For the t ^ dd final state we have considered all combinations of final state bottom 
and light down quarks. For the t —>• dl~^ final state, we have considered bottom quarks with 
all three lepton generations. For the t ^ uv final state, we studied top and charm quarks. 
For this final state, we have also considered the case where the LSP stop is lighter than the 
top in which case the stop decays to W~^hv via an off-shell top. 

The results are presented in Figs. 1-3. One can see that the whole region of natural stop 
masses < 800 GeV is excluded, except for purely hadronic prompt decay, t —?■ dd'. There 
is also some allowed parameter space by these searches for prompt decays with neutrinos, 
however we expect that prompt missing energy searches will further constrain some of the 
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Figure 2. 95% CL exclusion curves for stop decays to a bottom quark and a charged lepton. In 
the upper left plot we show t —>■ 6e+, in the upper right plot t —>■ and in the bottom 

plot we have t —)■ 6t+. While all displaced searches of ATLAS [28] and CMS [29] have been recast, 
we only display the curves yielding the strongest constraints. Stable searches for HSCPs [33] were 
recast for an unstable R-hadron. The prompt searches here correspond to leptoquark searches at 
the Tevatron and LHC [35-43]. 



Figure 3. 95% CL exclusion curves for stop decays to an up-type quark and a neutrino. In the left 
plot we have t ^ tv and in the right plot t —>■ cv. While all displaced searches of ATLAS [28] and 
CMS [29] have been recast, we only display the curves yielding the strongest constraints. Stable 
searches for HSCPs [33] were recast for an unstable R-hadron. The prompt searches here correspond 
to the R-parity conserving supersymmetric searches [3, 4, 44, 50, 51] 
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region. 

The bounds derived from the various combinations of bottom and light down quarks for 
the t ^ dd channel are extremely similar. One would expect the additional displacements 
of the bottom quarks to decrease the efficiency of vertex reconstruction with respect to the 
case of two light quarks. However, we found this effect is not highly signihcant (see Fig 
8 ). The reason is that the stop decay vertex contains tracks from the ii—hadron remnant 
along with radiation from the two bottom quarks before hadronization. These tracks can 
be reconstructed as a secondary vertex even if the H-meson decays are reconstructed as 
tertiary vertices. We expect this to be the case when a bottom quark is exchanged for a 
down quark for all stop and gluino channels in consideration. Thus, we have conhned our 
simulations for all other channels to include bottom quarks in the hnal states and have not 
considered the lighter d generations. For the Higgsino, which does not form an i?—hadron, 
the bounds for 3rd generation quarks are likely weaker than for light quarks which do not 
have secondary displacements, but the qualitative differences are expected to be small. 

The bounds derived from the t ^ tv and t ^ cv hnal states are also very similar 
except for the ATLAS DV+jets search. This difference is due to the fact that the top can 
decay with larger jet multiplicity than the lighter quarks, thus strengthening the bounds 
derived from hnal states with top quarks. However, since typically, the strongest bounds 
are derived from the ATLAS DV+MET search for hnal states with neutrinos, this is not of 
great signihcance. We expect the same to be true for the case of gluino and Higgsino LSPs 
and therefore consider only tops in the hnal states for their decays. 

The ATLAS DV+MET search bounds are presented only for the hnal states br, tv and 
cv, for which this search gives the strongest bounds (ATLAS DV+jets is comparable for 
the br hnal state). The bump in the bound at ~ 200 GeV LSP mass is a consequence of 
the MET > 180 GeV cut. 


4.2 Gluino LSP 

For the case of gluino LSP, we have considered the three types of hnal states as in Table 1 
corresponding to 5 —?■ tbb,ttv,tbl. Only direct production of the gluino has been simulated 
and the stop has been assumed to be decoupled for the gluino production cross section. 
The results are presented in Figs. 4-5. One can see that the whole region of natural gluino 
masses nig < 1500 GeV is excluded for long lived gluinos. There may be a small prompt 
region for the tbb hnal state that may still be somewhat natural. We expect the g —)• tbi 
hnal states to be constrained by same-sign dilepton and other relevant prompt searches 
which we have not attempted to recast here. 

As explained above, we have not considered the case of light down quarks in the hnal 
state, since the bounds for these are expected to be similar to those with bottom quarks. 
For the case of up type quarks in the hnal state, we only present results for top quarks, i.e. 
the channel g —?• ttv. As noted above, the decays g —)• tcv and g —>• tdv have very similar 
constraints for all channels. 
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Figure 4. 95% CL exclusion curves for gluino decays to tbb (left) and tiiy (right). While all 
displaced searches of ATLAS [28] and CMS [29] have been recast, we only display the curves yielding 
the strongest constraints. Stable searches for HSCPs [33] were recast for an unstable R-hadron. 
The prompt searches here correspond to the R-parity conserving supersymmetric searches [2, 5] 
and the RPV gluino search [45] . 



m-g (GeV) 

Figure 5. 95% CL exclusion curves for gluino decays to a top, a bottom quark and a charged 
lepton. In the upper left plot we have g —>■ tbr, in the upper right plot g —^ tbfi, while in 
the bottom plot we have g —)■ tbe. Although all displaced searches of ATLAS [28] and CMS [29] 
have been recast, we only display the curves yielding the strongest constraints. Stable searches for 
HSCPs [33] were recast for an unstable R-hadron. LHC searches for 3 or more leptons, same-sign 
leptons events, and/or multiple (6-)jets are not shown here, and are expected to constrain prompt 
decays. 
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Figure 6. 95% CL exclusion curves in the charged Higgino lifetime vs. mass plane, for Higgsino 
decays to a top, a bottom quark and a charged lepton. In the upper left plot we show H'^ —>■ tbb 
and —>■ bbb, in the upper right plot —>■ tbr and > hbr, and in the bottom plot we 

present > ttv and > btv. The ratio of neutral to charged Higgino lifetime depends on their 
mass and is indicated on all three plots. While all displaced searches of ATLAS [28] and CMS [29] 
have been recast, we only display the curves yielding the strongest constraints. Stable searches for 
HSCPs [33] were recast for an unstable Higgsino. LHC searches for 3 or more leptons, same-sign 
leptons events, and/or multiple (6-)jets are now shown above and are expected to constrain prompt 
decays. 


4.3 Higgsino LSP 

For the case of Higgsino LSP, we have considered direct prodnction of charged and nen- 
tral Higgsinos. We simnlated the prodnction of both particles simnltaneonsly nsing the 
corresponding mixing matrices for the relevant valne of /i ~ m^. Approximately 61% of 
prodnced Higgsinos are charged over the fnll range of masses, giving HSCP searches sen¬ 
sitivity to long-lived Higgsinos. The resnlts are plotted in Fig. 6 . We deconple all the 
snperpartners except for the stops, which we take to be degenerate. This prodnction mech¬ 
anism may be snbdominant to glnino or stop prodnction, which we have not considered. 
However, we obtain strong bonnds from direct prodnction alone, allowing for a conservative 
limit which rnles ont Higgsinos np to 800 GeV, except for a region of short lifetimes above 
~ 300 GeV. This region is expected to be constrained by prompt searches for large jet 
mnltiplicity events for [t/b)bb final states, same sign di-lepton for {t/b)bT final states and 
prompt MET searches for {t/b)tv final states as well as other prompt searches. Recasting 
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these searches is beyond the scope of this paper. 

We have simulated three channels, given in Table 1. For each coupling under consider¬ 
ation, the charged and neutral Higgsinos decay differently. For the r]" coupling, the decay is 
> tbb and —)• bbb. For the rj' coupling, the decay is either —)• tbr and —)• bbr 
or —)• ttv and —)• btu. In each case, the neutral (charged) Higgsino decays into 

an off-shell stop and a top (bottom) quark. The stop then decays via the RPV coupling. 
For down type quarks in the final state, we consider only bottoms quarks since the bounds 
for lighter generations are very similar. For up type quarks in the final state we consider 
only top quarks, again since the bounds for lighter generations are expected to be almost 
the same. Furthermore, we consider only r leptons in the final states since we expect the 
bounds for lighter generations to differ very slightly as was shown for the case of a gluino 
LSP. 

The charged and neutral Higgsinos have different lifetimes due to differences in phase 
space of the final states. The plotted lifetime is the lifetime. For below 200 GeV, 
the Higgsino decays via an off shell top. At 100 GeV, the neutral Higgsino is effectively 
stable due to phase space suppression. The qualitatively different features of these results 
are due to the presence of two different average lifetimes in each data sample. For example, 
in Fig. 6, the ATLAS DV+MET search becomes important for masses where the neutral 
Higgsino becomes long-lived and contributes to MET. Prompt lifetimes are ruled out at 200 
GeV due to the presence of longer lived neutral Higgsino which decay within the tracker. 

5 Conclusions 

We have considered the experimental bounds on supersymmetric theories with R-parity vi¬ 
olation and a long-lived LSP. The main searches providing the constraints are the GMS 
displaced dijet search [29] and the ATLAS search for multitrack displaced vertices + 
/x/e/jets/MET [28], which we have fully recast and applied all cuts and vertex recon¬ 
struction procedures. In addition we also considered the GMS HSGP search [33], as well 
as some of the leading prompt searches whose exclusion bounds can be directly applied 
without recasting. Our main results are the exclusion plots in the misp ~ plsp plane in 
Figs. 1-6, where the LSP is a stop in Figs. 1-3, a gluino in Figs. 4-5, or a Higgsino in Fig. 6. 
The various plots for a given LSP correspond to different decay channels either from the 
ordinary holomorphic RPV operator of (1.1) or the non-holomorphic dRPV operators of 
(1.2). All regions with significant displacement, cr > 1 mm, are excluded for the natural 
SUSY mass regions < 800 GeV, rUg < 1500 GeV or < 800 GeV. The only significant 
unconstrained regions are those corresponding to prompt stop decays to dijets or prompt 
Higgsino decays. 
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Appendix 

A Data and Simulations 
A.l Monte Carlo Tools 

Parton level events are generated using the MadgraphB event-generator [52] , utilizing model 
files created by Feynrules [53]. Hadronization, parton showering, and i?-hadron physics 
are simulated using Pythia 8 [54], which allows the eap-y color flow associated with the qqd* 
and udd operators. Long-lived particles are displaced before passing events from Madgraph 
to Pythia. For the detector simulation, we use Delphes 3 [55] and ROOT 6.02 [56] with 
the default ATLAS and CMS detector geometries, resolutions, and efficiencies (excluding 
tracking efficiency). The Delphes 3 detector has an inner detector (ID), electromagnetic 
and hadron calorimeters (ECAL/HCAL), and a muon system (MS). Charged SM particles 
are propagated in a uniform magnetic field in the tracker. FastJet [57] is used to cluster 
jets according to track and calorimeter information from the Delphes particle-flow recon¬ 
struction. Tracking vertex resolution is taken to be perfect, but nearby vertices are merged 
in the vertex reconstruction procedures. 

A.2 Tracking Efficiency 

The standard tracking simulation provided in Delphes 3, which assumes tracking efficien¬ 
cies for particles originating from the PV, provides a poor reproduction of the actual tracking 
efficiency for tracks originating from displaced vertices. The track reconstruction efficien¬ 
cies in ATLAS [58] and CMS [59, 60] drop off rapidly for vertices displaced from the origin. 
For the ATLAS and CMS displaced vertex searches presented in this work, offline tracking 
techniques (“re-tracking”) are used to regain some of the lost efficiency. Even with displaced 
tracking algorithms, the probability to reconstruct a track in the CMS or ATLAS inner de¬ 
tectors depends on the track’s orientation and the displacement of its production vertex. 
This effect is incorporated in our analysis by assigning displaced tracks efficiencies based on 
simulated displaced track reconstruction efficiencies of the ATLAS and CMS detectors. For 
the ATLAS detector, we use an efficiency parametrized by the transverse impact parameter, 
do, and transverse momentum, px, of each track. In addition to tracking efficiency falloff 
with do, ATLAS vertex reconstruction efficiency also drops off steeply at a radial distance 
corresponding to the second pixel barrel layers of the ID. To account for this effect, we add 
an additional factor of 75% to the track reconstruction efficiency for tracks which originate 
outside of the second pixel layer. The CMS track reconstruction efficiencies are simulated 
under Run 1 tracking conditions [59] and are parametrized as a function of production 
radius. A cut on track transverse impact parameter do < 30 cm is imposed, motivated by 
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displaced tracking efficiencies for cosmic muons [60]. The final reconstruction efficiencies 
do not strongly depend on the the shape and scale of the tracking efficiency nor on their 
parameterization either as a function of production radius or impact parameter. 

Tracks of charged i?-hadrons which decay in the tracker are assumed to not recon¬ 
struct and are neglected in this study. It is possible that some of these tracks could be 
reconstructed by ATLAS or CMS, but this depends on the tracker hits registered in the de¬ 
tector before the charged i?-hadron decays and the details of the tracking algorithm. Tracks 
of charged i?-hadrons which decay in the beam pipe do not register hits in the detector 
and therefore are not reconstructed. i?-hadrons decaying in the tracker may have tracks 
which exhibit multi-prong decays or which disappear altogether, and may also not recon¬ 
struct properly depending on the reconstruction algorithm. The track from the charged 
i?—hadron would appear to be prompt because it points back to the primary vertex (PV). 
The ATLAS multitrack DV search ignores prompt tracks with transverse impact parameter 
do < 2 mm in reconstructing displaced vertices. The CMS dijet search allows for only one 
prompt track to be associated to each displaced jet, so the search could retain sensitivity to 
events in which the charged i?-hadron tracks are reconstructed. Therefore, it is reasonable 
to neglect the track of the i?-hadrons in this study. The curvature of charged i?-hadron 
trajectories, typically small due to the large momentum, is neglected. 

B Searches 

B.l ATLAS DV+/i/e/jets/MET 

In this section, we describe the ATLAS [28] multitrack search for displaced vertices along 
with our procedure. Trigger and selection cuts vary between the different final state topolo¬ 
gies while displaced vertex reconstruction remains the same for all cases. 

B.1.1. Trigger Requirements 

For the DVTmuon search, a trigger muon is required to have: 

• Pt > 50 GeV 

• pseudo-rapidity \rj\ < 1.07 to be within the MS acceptance 

• tracker hits both in the ID and MS 

For the DVTelectron search, a photon trigger is implemented since electrons with large 
transverse impact parameters may not have a measured track in the inner detector and 
may be identified as photons. The search triggers on either one or two photons with large 
transverse momentum. The requirement is: 

• either one photon with pT > 120 GeV or two photons, each with pT > 40 GeV 

For the DVTjets search, jets are reconstructed using the anti- A:* algorithm with size pa¬ 
rameter Ai? = 0.6. The trigger requires: 

• four jets with pT > 80 GeV or five jets with px > 55 GeV or six jets with px > 45 GeV 
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The trigger requirement for the DV+MET search is: 

• > 80 GeV 

B.1.2. Event Selection 

Offline cuts filter displaced tracks which undergo a re-tracking procedure prior to hnal 
analysis. The details of these cuts and the re-tracking are given in [28]. Following this 
procedure, hnal event selection cuts are implemented for the various channels. In our 
analysis, we require the hnal event selection criteria on the reconstructed leptons, jets, and 
MET since they are more restrictive than the trigger requirements in all cases. The trigger 
leptons must have a truth level origin within the ID. 

All events must have a reconstructed DV with track multiplicity > 5 and invariant 
mass niDv > 10 GeV, where each track is assigned a charged pion mass for the invariant 
mass determination. 

For the DVTmuon search, the cuts are similar to the muon trigger, requiring: 

• Pt > 55 GeV 

• |??| < 1.07 

The DV+electron search cuts require the electron candidate to have: 

• Pt > 125 GeV 

• \r]\ < 2.48 

Furthermore, for both lepton searches, the lepton track must be displaced, corresponding to 
a transverse impact parameter which satisfies do > 1.5 mm and is required to pass within 
0.5 mm of a reconstructed DV. 

Guts for the DV+jets search are stricter than the trigger requirements requiring: 

• four jets withpr > 90 GeV or five jets withpT’ > 65 GeV or six jets withpy > 55 GeV 
For the DV+MET search, the requirement is: 

• > 180 GeV 

There is a subtlety regarding whether LSPs which escape the tracker contribute to the 
transverse energy of the event for the DV+MET search. A charged i?-hadron or Higgsino 
which decays outside of the entire detector may be reconstructed as a muon from its tracks 
in the MS and therefore be included in the energy of the event. LSPs, charged or neutral, 
which decay outside of the ID but within the calorimeter or within the MS may contribute 
to the energy of the event either from calorimeter deposits or from the tracks of its decay 
products. In our analysis, we treat such LSPs as follows: Gharged LSPs which traverse the 
entire detector and leave a track in the ID are reconstructed as muons and the transverse 
momentum of the LSP’s track does not contribute to MET. Transverse energy deposited by 
the decay products of either charged or neutral LSPs which decay within the calorimeter 
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is removed from the MET of the event with the exception of neutrinos in the final state. 
For LSPs decaying in the MS, charged decay products are not included in the MET of 
the event. We neglect the calorimeter deposits from the i?-hadrons themselves since they 
are expected to be small. This MET reconstruction procedure for decays in the various 
regions of the detector is a conservative approach since, in reality, some of this energy is 
not reconstructed and will contribute to the MET of these events. 

B.1.3. Vertex Reconstruction 

Events satisfying the trigger requirements must also contain at least one reconstructed 
DV. Vertex reconstruction in the ATLAS search is performed considering only displaced 
tracks with do > 2 mm and pT > 1 GeV. All tracks are extrapolated in the direction 
opposite to their momentum. Tracks which pass the criteria are used to reconstruct DVs 
by means of an algorithm detailed in [28] which clusters seed vertices iteratively. A seed 
vertex is defined as two coinciding tracks which each satisfy the criteria d ■ p > —20 mm, 
which ensures consistency between the position of the vertex and momentum of the track. 
Here d is defined as the distance vector between the position of the DV and that of the first 
primary vertex and p is the seed track’s momentum. After the tracks are clustered into 
vertices, any reconstructed DVs within 1 mm are merged into a single vertex. 

After vertex reconstruction, DVs are further required to have transverse distance Lxy < 
300 mm and \zdv\ < 300 mm with respect to the PV. In order to minimize background 
from tracks originating at a PV, each DV must also have a transverse distance of at least 
4 mm from any primary vertex. Finally, DVs situated within dense regions of the detector 
are vetoed using a 3D map of the transverse plane of the detector [28]. 

For our analysis, vertex reconstruction is identical for all the final state searches. Tracks 
considered for vertex reconstruction are required to have pT > 1 GeV and impact parameter 
do > 2. The truth level origins of all tracks must also satisfy 4 mm < Lxy < 300 mm and 
|z| < 300 mm. Our vertex reconstruction algorithm searches through all tracks that satisfy 
these criteria and iteratively clusters their origins. First, two track origins which are within 
a 1 mm radius from one another are grouped together as a DV. The DV position is defined 
as the average position of all of the track origins in the group. DVs are then combined if 
their distances are less than 1 mm apart. Furthermore, at least two tracks in each DV are 
required to satisfy d ■ p > —20 mm in consistency with the seed vertex track requirement 
detailed above. DVs are vetoed at positions in the transverse plane which are mapped by 
the ATLAS study [28] as dense material regions of the detector. 

For the DV+lepton channels, lepton tracks are extrapolated in the opposite direction 
of their momenta. The extrapolated track is then required to have a maximum distance 
of closest approach to the reconstructed DV of 0.5 mm. This procedure has important 
implications for signals with final state bottom or charm quarks since these decays result in 
additional displacement from the LSP decay vertex, and the DV clustering algorithm could, 
in principle, reconstruct these vertices separately. The ATLAS study associates each track 
with only one DV. However, the lepton track is not required to fit to the reconstructed 
DV and may therefore pass within 0.5 mm from multiple vertices, as in the case of the 
decay products of very boosted bottom or charm quarks. This is taken into account in our 
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Figure 7. Example event-level efficiencies for the ATLAS displaced vertex searches. Efficiencies 
are shown for = 200 (red), 300 (green), 500 (orange), 1000 (blue) GeV. Efficiencies generally 
increase with increases stop mass. 

procedure since truth level tracks are always associated with only one DV and lepton tracks 
are extrapolated and not required to be directly associated with the DV. 

B.1.4- Efficiencies 

With these assumptions for vertexing and tracking efficiency, we reproduce the life¬ 
time dependence of the efficiencies for the models and channels provided in [28]. For the 
DV+jets and DV+MET channels, only event-level efficiencies are calculated since the 
cuts for these channels are required at the event level and not at the vertex level. For 
the DV+lepton channels, vertex-level efficiencies are calculated and then translated into 
event-level efficiencies. We accept up to two DVs per event and calculate a DV level effi¬ 
ciency, eov- III consistency with the ATLAS study, for the DV+lepton channels, the event 
level efficiency, Cev, is defined as the probability for a signal event containing two DVs to 
satisfy all selection criteria for at least one of the DVs in the event. This definition relates 
to the vertex level efficiency, through, 

^ev = 2e£)y ~ ^DV- (I^-I) 

Example event level efficiencies are shown in Fig. 7 and have been used in our calculation 
of the exclusion plots in Sec. 4. 
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B.2 CMS Displaced Dijet 


CMS [29] performs two searches for signals with high and low average transverse displace¬ 
ments. We recast the study using the high-(La;y) selection criteria, which is strong enough 
to place limits on the full range of displaced lifetimes for the models considered in this 
paper. 

B.2.1. Trigger Requirements and Offline Cuts 

CMS uses a dedicated displaced-jet trigger for this search which requires Ht > 300 GeV, 
where Ht is the sum of the transverse energy of all the jets in the event with pT > 40 GeV 
and \r]\ < 3, and at least two displaced jet candidates (isolated leptons are also considered) 
each satisfying: 

• Pt > 60 GeV and \r]\ < 2, 

• at most 2 associated tracks with three-dimensional impact parameter less than 300 pm, 

• at most 15% of the jet energy is carried by associated tracks with transverse impact 
parameters less than 500 pm. 

At trigger level, jets are clustered according to calorimeter information only. 

The offline event selection cuts are similar to the triggers. The total transverse energy 
requirement is raised to Ht > 325 GeV for the fully reconstructed jets, using anti-ZcT algo¬ 
rithm with AR = 0.5. For our jet reconstruction, we use track and calorimeter information 
from the Delphes particle-flow reconstruction, and cuts are imposed according to the offline 
selection criteria made by GMS. 

B.2.2. Dijet Reconstruction and Final Event Selection 

Every jet pair is checked for consistency with the displaced dijet hypothesis. Displaced 
tracks, defined by transverse impact parameter do > 0.5 mm, associated to each jet pair are 
grouped together and fitted to a common secondary vertex. Tracks originating at a radial 
distance greater than 50 cm have zero acceptance because the strip tracker layers no longer 
contain 3D information, and the search requires track seeds to be within the part of the 
tracker with stereo tracking information. The full analysis uses an adaptive vertex fitter 
which is not easily implemented. We mimic the vertex fitter by merging the origin of two 
tracks if their truth-level origins are within a distance of 1 mm. If more than one displaced 
vertex is found between the two jets, we take the one with largest track multiplicity. Each 
of the two jets must have at least one track originating from the secondary vertex, which 
must be significantly displaced from the primary vertex of the event. The search considers 
the transverse displacement of a secondary vertex, Lxy, to be significant if Lxy > 8 (TL^cyi 
where ctl^cv the uncertainty in Lxy The primary vertex resolution for GMS is 12 microns 
in each dimension [61], which is negligible compared to the uncertainty in the secondary 
vertex. We take the uncertainty in Lxy be a constant 300 pm for all DVs as a conservative 
choice, requiring the displaced vertex to be at a distance of 2.4 mm away from the primary 
vertex in the transverse plane. For decays with short decay lengths (for which this cut is 
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relevant), the resolution of L^y is better than 300 /rm given the good resolution near the 
center of the detector. This cut does not signihcantly affect the hnal efficiencies because 
decays which occur this close to the PV begin to have tracks which point back to the PV 
and fail the requirements of the displaced-jet trigger. 

For hnal event selection, further selection criteria are imposed for each displaced dijet 
candidate passing the above criteria. Clusters of maximal track multiplicity are formed 
from displaced tracks associated to dijet candidates based on the transverse displacements 
of each individual track, is determined by hnding the intersection of the dijet 

transverse momentum (the sum of the two reconstructed jet pt) with the trajectories of 
the individual tracks. In the full analysis, the individual tracks are clustered using 

an algorithm with size parameter 0.15 Lxy For this recast, we approximate this procedure 
by grouping tracks with within a transverse distance 0.15 Lxy together in a cluster. 

The associated track cluster must also satisfy: 

• at least one track from each jet belongs to the cluster 

• invariant mass of cluster > 4 GeV 

• sum of track py > 8 GeV 

Tracks are assigned the mass of the charged pion for the purpose of determining the invariant 
mass. 

Finally, a background discriminant is formed by the secondary vertex track multiplic¬ 
ity, the cluster track multiplicity, the cluster RMS of and the fraction of secondary 

vertex tracks having positive signed impact parameter. We neglect the RMS and signed 
impact parameter discriminants since we do not implement the clustering algorithm used 
by the search and the vertex and cluster track multiplicity affects the value of the full dis¬ 
criminant most signihcantly because background DVs typically have low track multiplicities 
compared to signal DVs. We use the distributions given in [29] to form the vertex and clus¬ 
ter multiplicity discriminant. The hnal event selection for the high (Lxy) search requires 
the background discriminant to be less than 0.8, prompt track energy fraction for each jet 
less than 9%, and no more than one prompt track associated to each jet. 

B.2.3. Efficiencies 

The high {Lxy) search observed 1 data event with an expected background of 1.13 ± 
0.15 ± 0.50. This places the excluded number of signal events for 95% conhdence level at 
3.7. The results are interpreted in the context of two models: (1) a heavy scalar particle, 
Hq, decaying into a pair of long-lived neutral Xq particles, which then decay into dijets, 
and (2) squark pair production, in which each squark decays via q —>• g(x5 udp). The 
study considers Xq and Xj* proper lifetimes from 0.1 to 200 cm. We match the reported 
efficiencies results within 20% over the full range of lifetimes for both 2-body and 3-body 
hnal states. 

Example event level efficiencies are shown in Fig. 8 and have been used in our calculation 
of the exclusion plots in Sec. 4. 
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Figure 8 . Example event-level efficiencies for the CMS displaced dijet search. Efficiencies are shown 
for rrij = 200 (red), 300 (green), 500 (orange), 1000 (blue) GeV. Efficiencies generally increase with 
increases stop mass. Due to the additional displacement of the 6-jets, the efficiencies are lower in 
the t ^ bd decay channel than in decays to d- and s-quarks. 

B.3 CMS Heavy Stable Charged Particles 

If the LSP lifetime is sufficiently long, searches for heavy stable charged particles (HSCP) 
may be sensitive to the scenarios considered in this paper. Here we briefly describe the 
CMS search presented in [33]. 

Triggered events must have either a reconstructed muon with px > 40 GeV measured 
in the ID or a large missing transverse energy > 150 GeV. The LI muon trigger can 

accept slowly moving particles arriving in the MS within the 25 ns of the proton bunch 
crossing or within the following 25 ns time window before the next bunch crossing. Neutral 
i?-hadrons or HSCPs that become neutral do not leave tracks in the muon system. Conse¬ 
quently, they are rejected by the online particle flow algorithm and are not reconstructed 
as objects. However, since it only leaves 10-20 GeV of energy in the calorimeter, there is 
large large in the event. In such cases, the HSCPs are still identified by the anoma¬ 

lously high energy loss in the inner tracker. The offline data selection requires a track with 
Pt > 45 GeV, |7y| < 2.1, and a dE/dx > 3 MeV/cm. 

It is straightforward to reinterpret the search results for unstable i?-hadrons. The 
excluded cross sections derived from the full CMS analysis already includes the acceptance 
and reconstruction efficiencies for stable i?-Hadrons. The upper limits for stop and gluino 
production cross section are rescaled by the percentage of events passing the pT and \r]\ cut 
that survive the length of the detector. The cut on dE/dx is neglected due to our limited 
detector simulation. Efficiencies for the colored particle to decay outside the detector is 
found by simulating parton-level events for stop production with masses between 100 GeV - 
1 TeV and gluino production with masses between 100 GeV - 1.5 TeV. The detector is taken 
to be a cylinder with radius of 7m and 11m half-length [62], and the long-lived particle is 
required to survive the entire detector distance. This is a conservative choice as it neglects 
stops that decay in the muon tracker which may still be picked up by the full analysis due 
to the presence of high pT tracks in the MS. 

The results are interpreted in terms of stable gluinos and stops, placing upper limits 
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of 1322 (1233) GeV and 935 (818) GeV in the cloud (charge-suppressed) model, respec¬ 
tively. GMS performs a “tracker-only,” “tracker+time-of-flight,” and “muon-only” search. 
The “tracker-only” search is sensitive to standard i2-hadrons and the charge-suppressed 
scenario because it does not rely on a charged i?-hadron tracks in the muon system, so we 
choose to recast this search for the case of unstable i?-hadrons. For charged higgsinos, we 
apply the “tracker+time-of-flight” search. 

B.4 Prompt Searches 

For our analysis, prompt searches have not been recast. Instead we have studied a partial 
subset of prompt searches whose bounds can be directly applied to the scenarios of interest. 
Recasting other searches may possibly result with more competitive limits than the ones 
considered here: 

• t —>• dd. The main prompt searches are the paired dijet searches performed at 
GMS [34]. Gonstraints for light quark searches are taken to be valid for 6-quark 
final states. 

• t —>■ d^'^. We consider leptoquark searches from the Tevatron and the LHG. Searches 

for leptoquarks are performed for "generation" leptoquarks, i.e., searches are per¬ 
formed for de^ [35-37], [38-40], and hT~^ [41-43] decays of the leptoquark. Gon¬ 

straints for light quark searches are taken to be valid for 6-quark final states. 

• t —)• ci>. We consider ATLAS supersymmetry search for charm squarks, c —>■ cN, 
where N is the LSP [44]. The applicable bounds are those for = 0 GeV. These 
bounds are much stronger than the monojet search [63]. 

• 5 —)• ttu. We study LHG supersymmetry searches for g —>• ttN. The applicable 
bounds are those for = 0 GeV. The strongest bounds are derived from the 0-1 
leptons and (> 3) 6-jets search from ATLAS [5] and 1 lepton and (> 6) jets search 
from GMS [2]. 

• i ^ W. We consider LHG supersymmetry searches for i —)• tN, where N is the 
LSP. The applicable bounds are those for = 0 GeV. The strongest bounds are 
derived from the 1-lepton searches from GMS [51] and ATLAS [50]. For < mt the 
bounds are taken from [3] and for mt < m^ < 200 GeV the bounds are taken from 
the measurement of spin correlations in ti events [4]. 

• 5 —)• tiv. We study LHG supersymmetry searches for g —>■ tiN. The applicable 
bounds are those for = 0 GeV. The strongest bounds are derived from the 0-1 
leptons and (> 3) 6-jets search from ATLAS [5] and 1 lepton and (> 6) jets search 
from GMS [2]. 

• 5—7- tbb. We study the ATLAS RPV search for g —?■ tbs [45], and take the search to 
be valid for additional 6-jet multiplicity. 
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Prompt searches may also apply to LSPs with large enough boosts such that its decay 
products point back to the PV. Specific displaced tracking algorithms must be used to 
reconstruct tracks with large transverse impact parameter. Without simulating the full 
detector and tracking algorithms, it is not known how the prompt search bounds extend 
for longer lifetimes and we therefore do not consider this possibility. 
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